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We investigate the finite-field ground state of tlie 5 = 1 antiferromagnetic-ferromagnetic 
bond-alternating chain described by the Hamiltonian Ti. — '}2e{S2e-i ■ S2i-'rJS2i ■ 821+1] + 
D {Si)'^—H Si, where J < and — cxd < D < 00. We find that two kinds of magnetization 
plateaux at a half of the saturation magnetization, the ^-plateaux, appear in the ground-state 
magnetization curve; one of them is of the Haldane type and the other is of the large- D-type. 
We determine the ^-plateau phase diagram on the D versus J plane, applying the twisted- 
boundary-condition level spectroscopy methods developed by Kitazawa and Nomura. We also 
calculate the ground-state magnetization curves and the magnetization phase diagrams by 
means of the density-matrix renormalization-group method. 

KEYWORDS: S=l antiferromagnetic-ferromagnetic bond-alternating chain, numerical methods, ground- 
state magnetization curve, i-plateau phase diagram, magnetization phase diagram 
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1. Introduction 

There has been a considerable current interest in the 
study of the ground-state properties of various one- 
dimensional quantum spin systems in a finite magnetic 
field. In this paper we discuss the case of the 5=1 bond- 
alternating chain, for which we express the Hamiltonian 
in the following form: 

~ Ti-cx + Hz , (1) 
Wcx = ^^{'^2«-i • + JS2e ■ S2e+i^ 

nz = -Hj2^S^,, (3) 

where Si is the S—1 operator at the £th site; J 
(— cx)< J<cxd) is the parameter representing the bond 
alternation of the nearest-neighbor interactions; D 
(— oo<Z)<oo) is the uniaxial single-ion-type anisotropy 
constant; H {H>0) is the magnitude of the external 
magnetic field applied along the z-direction. Hereafter, 
we denote by M the z component of the total spin 
5tot = , where N, being assumed to be a mul- 

tiple of four, is the total number of spins in the system. 

The finite-field as well as the zero-field ground-state 
properties of the present system in the case of the 
antiferromagnetic-antiferromagnetic bond-alternating 
chain with J > have already been investigated by 
using mainly numerical methods.^' ^ As for the zero- field 
properties, the ground-state phase diagram on the J 
versus D plane, in which the Neel, Haldane, large-!?, 
and dimer phases appear, has been determined.^' ^ On 
the other hand, the results for the ground-state magne- 
tization curve show that the magnetization plateau at a 
half of the saturation magnetization, which is called the 
i-plateau, appears for arbitrary values of J except for 
J=l at least when D>Q.^ 



In this study we explore the finite-field ground- 
state properties of the system in the case of the 
antiferromagnetic-ferromagnetic bond-alternating chain 
with J<0, focusing our attention mainly upon the ^- 
plateau appearing in the magnetization curve. When 
D = 0, the present system is mapped onto the isotropic 
S~2 uniform antiferromagnetic chain in the limit of 
J —>-~oo. Thus, it is considered that, no plateau appears 
in the ground-state magnetization curve in this limit 
when D — 0. As we have already discussed,^ on the other 
hand, the ^-plateau appears when J — and D>—^. 
These imply that the '^-plateau'-'no-plateau' transition 
occurs at a finite value of J for a certain region of D in- 
cluding D = 0. Furthermore, as has already been shown, 
the magnetization plateau at one third of the saturation 
magnetization, the ^-plateau, appearing in the ground- 
state magnetization curve of the 5* = | uniform antiferro- 
magnetic chain with single-ion-type anisotropy changes 
its character from the Haldane-type to the large- D-type 
at a finite value of the anisotropy constant. We expect 
that this type of 'Haldane-type-^-plateau'-'large-D-type- 
^-plateau' transition may take place also in the present 
5 = 1 antiferromagnetic-fcrromagnetic bond-alternating 
chain. To discuss quantitatively these transitions and to 
complete the ^ -plateau phase diagram on the D versus 
J plane are our main purposes. We also aim at calcu- 
lating the magnetization phase diagram as well as the 
ground-state magnetization curve. 

In the next section (§2) we discuss the ^-plateau phase 
diagram. In determining the phase diagram, we apply the 
twisted-boundary-condition level spectroscopy (TBCLS) 
analyses of the numerical diagonalization data, devel- 
oped by Kitazawa'* and also by Nomura and Kitazawa.^ 
Then, we discuss in §3 the ground-state magnetization 
curve, which we calculate by using the density-matrix 
renormalization-group (DMRG) method proposed origi- 
nally by White. ^ We also discuss the magnetization phase 
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Fig. 1. Schematical pictures representing (a) the 'Haldane-type- 
j-platcau' phase and (b) the 'large- D-type-^-plateau' phase. 
The thick solid and dotted lines stand, respectively, for the an- 
tiferromagnctic and ferromagnetic bonds, and the small open 
circles for the S=^ spins. Each large open circle surrounding 
two S = 5 spins stands for an operation of constructing an S=l 
spin by symmetrization. The 15 = ^ spin with an upward arrow 
is in the up state f- Furthermore, two S=^ spins in a flat el- 
lipse form a singlet dimer (Ti ^ it) > E^nd two S=^ spins in a 
rectangular form a triplet dimer --i=(l"J.-|-J.|). 



diagram which is obtained from the resiih for the mag- 
netization curve. Finally, concluding remarks are given 
in §4. 

2. ^-Plateau Phase Diagram 

Prom the considerations presented in the previous 

section, it is expected that three phases, i.e., the 'no- 
plateau', 'Haldane-type-i-plateau', and 'large- D-type- 
|-plateau' phases appear in the ^-plateau phase dia- 
gram. Schematical pictures in terms of a kind of the 
valence-bond-solid picture, which represents the latter 
two phases, are shown in Fig. 1. 

Wc can determine the transition lines between two 
of the above thee phases rather precisely by applying 
the TBCLS methods in the following way. Both the 
'Haldane-typc-i-platcau'-'no-platcau' transition and the 
'large- D-type-i-plateau'-'no-plateau' transition are the 
Berezinskii-Kosterlitz-Thouless transition^ accompany- 
ing no spontaneous breaking of the translational symme- 
try of the Hamiltonian H. Therefore, the critical points 
j^H,no) j^LD.no) ^£ ^^jch are, respectively, for the 
former and latter transitions, for a given value of D (or, 
Dc^'"°^ and of D for a given value of J) can be 

estimated by means of Nomura and Kitazawa's TBCLS 
method. We denote by -Bq {N, M) the lowest energy 
eigenvalue of Tiox with the periodic boundary condition 
within the subspace determined by the values of N and 
M, and also by ^^q"^' {N, M, P) the lowest energy eigen- 
value of Ticx with the twisted boundary condition within 
the subspace determined by the values of N, M, and P, 
where P is the eigenvalue of the space inversion opera- 
tor. Then, Nomura and Kitazawa's TBCLS method^ pre- 
dicts that Jc^'""-* (£)c^'"°') is obtained by extrapolating 



Jc^''^°\n) [D^''^°\n)) determined from the equation 
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<)(iV,f + 2)+Er(A^,f-2)] (4) 



to N^oo. Similarly, J^'"^'""^ (^D^^'^°^) is obtained 
by extrapolating Jc^^'"°^ (AT) (£)c'"°'"°^(-^)) determined 
from 
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(5) 



to A^^oo. 

On the other hand, the 'Haldane-type-|-plateau'- 

'largc-_D-typc-i-platcau' transition is expected to be of 
the Gaussian type as in the case of the 'Haldane-type- 
i-plateau'-'large-f-type-l-plateau' transition in the the 
S* = I uniform antifcrromagnctic chain with single-ion- 
type anisotropy,"^ discussed in the previous section. Thus, 
we can apply Kitazawa's TBCLS method^ to estimate 
the critical point J^^'^^^ of J for a given D (or, oi^'^^^ 
for a given J) in this transition. According to this 
method, Ji«'"°) (i?(«-"°') can be estimated by extrap- 
olating ji"'^°)(Af) {Di^'^^\N)) determined from 



E, 



(T) 



E, 



(T) 



(6) 



to N^oo. It is noted that in the limit of N—>-oo, 



the value of E, 



(T) 




A^, -J, — 1 ) is smaller or larger than 



that of E^^ (^N, -^,-1-1^ depending upon whether in the 

'Haldane-type-i-plateau' region or in the 'large- D-type- 
^-plateau' region. 

Determining the ^-plateau phase diagram on the D 
versus J plane, we have first calculated numerically 
the energy eigenvalues E^ (^N, f +2) , E^^^ (^N, f - 2) , 

EP(^N,f,+iy and 4"^^ (aT, f , -l) for a vari- 
ety values of J and D for finite-size systems with 
N = 8, 12, and 16 spins, employing the computer 
program package KOBEPACK^ coded by means of 
the Lanczos technique. Then, we have solved numer- 
ically eqs. (4), (5), and (6), respectively, to evaluate 



(H,no) 



{N) (D^"'"°^(Af)), J, 



(LD,no) 



(LD,no) 



and j!:^'^^\n) {D^J^'^^\n)) . Finally, we have ex- 
trapolated these finite-size values to N~^(x to esti- 
mate (i^(«^"°^ (d(^°'"°^ and ./i^'^^' 

(D^"'^°^) by assuming that the A"-dcpendcnccs of 
the finite-size values are quadratic functions of N~'^. 
Plotting the extrapolated results on the D versus J 
plane, we have obtained the ^-plateau phase diagram 
as depicted in Fig. 2. It is noted that, for example, 
j^(H,no) Qg _|_ Q Qi fQj. £) = o. Furthermore, the three 

transition lines in the phase diagram meet at the point 
(J= -8.13 ± 0.01, £> = 1.19 ± 0.01). 

As can be seen from Fig. 2, the 'large- D-typc-^- 
plateau'-'no-plateau' transition line is almost parallel to 
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^-plateau phase diagram on the D versus J plane. 

NO, H, and LD stand, respectively, for the 'no-plateau', 

1 



Fig. 2. 
Here, 

'Haldane-type-^-platcau', and 'large- D-typc-^-platcau' regions. 
The three transition lines meet at the point (J = —8.13 ± 0.01, 
D = 1.19 ± 0.01), which is marked by the double circle. 



the J-axis in the whole region of J<— 8.13. This re- 
sult may be attributed to the fact that the present 
5 = 1 system in this region is almost equivalent to the 
5 = 2 uniform antiferromagnctic chain with single-ion- 
type anisotropy, since | J| is fairly large with J<0. Thus, 
we may expect a direct transition from the 'no-plateau' 
phase to the 'large- D-type-|-plateau' one in the above 

i -plateau in 



5 = 2 chain, in contrast to the case of the 
the corresponding 5 = | chain,'' mentioned in §1. Fig- 
ure 2 also shows that the 'large-D-type-^-plateau' phase 
appears when both | J| ( J < 0) and D (D > 0) are suffi- 



ciently large. The fitting of several values of jj;^'^^^ 

sufficiently large D's to a quadratic function of D^^ 

1 



for 



sug- 
gests that the 'Haldane-type-^-plateau'-'large-L'-type-^- 
plateau' transition line approaches J = —2.05 ± 0.05 in 
the limit of D— >oo. 

3. Magnetization Curve and Magnetization 
Phase Diagram 

Let us begin with discussing the saturation field Hg for 
the present system in the case of J<0. It is given by 

H, = {N, N)-E^^\n,N-1) = D + 2 (7) 



W AT) - eP{N,N- 2)} = ^E^ 



or 

Hs = l- Km 

(8) 

depending on whether 2(_D + 2) > Eb or 2{D + 2) < Eb, 
where Eb is the maximum energy of the ferromagnetic- 
two-magnon bound state. In principle, the energy Eb 
can be calculated analytically. However, this calculation 
in the case of J 7^ 1 has not yet been done as far as we 
know,^" and therefore we estimate Hs numerically when 
2{D+2)<Eb. 

We have applied the DMRG method^ to calculate the 
ground-state magnetization curve, which we define here 
as the average magnetization m(N)(=M/N) per spin 
versus the reduced field H/Hs curve. In this application, 
we usually have to impose open boundary conditions 
for the Hamiltonian treated. In the present system, we 
have two kinds of open boundary conditions, Oa and Of, 




^(Ox), 



Fig. 3. Ground-state magnetization curves in the N^oo limit 
obtained for the cases of (a) J=— 1.0 and D = 3.0 and of (b) 
J = -10.0 and 13 = 3.0 



which are, respectively, the cases where the edge bonds 
are the antiferromagnctic and ferromagnetic bonds. We 
adopt, respectively, the Oa and Of conditions when we 
calculate the magnetization curves in the 'Haldane-type- 
^-plateau' and 'large-D-type-^-plateau' regions. This is 
because, as can be seen from the schematical pictures 
for these phases shown in Fig. 1, no degrees of freedom 
appear at both edges of the chain in the above adoption, 
which leads to a less boundary effect. 

We have calculated the lowest energy eigenvalue 
^ (TV, M) (x=f or a) under the open boundary con- 
dition Ox of Hex within the subspace determined by the 
values of N and M for AT = 36, 64, and 84 and for all 
values of M (M = 0, 1, • • • , N). From these results the 
ground-state magnetization curves for TV = 36, 64, and 
84 can be readily obtained.^ The resulting magnetiza- 
tion curve for each TV is a stepwisely increasing function 
of H, starting from m{N) = Q at H = Ho{N) and reach- 
ing to the saturation value ■m{N) = l at H = Hs{N). We 
note that iJg = limjv^oc Hs{N). 

Except for plateau regions, a satisfactorily good ap- 
proximation to the magnetization curve in the N^oo 
limit, i.e., the ?7?,(= 777(00)) versus H/Hg curve, may 
be obtained by drawing a smooth curve through the 
midpoints of the steps in the finite-size magnetization 
curves.^ As for the estimations of 7Jo(=-ffo(oo)), which 
is nothing but the energy gap between the ground state 
and a first excited state in the case of H = 0, as well 
as of Hi and H2, which are the lowest and highest val- 
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Fig. 4. Magnetization phase diagrams on the H/Hs versus J 
plane obtained for the cases of (a) £' = 0.0 and (b) I? = 3.0, where 
Hq I H (open squares) , H\/H (soHd circles) , and H2 / H (open cir- 
cles) ate plotted as functions of ,7. In the regions A, Bx (X=H 
or LD), and C, the average magnetization m per spin is equal 
to (the 0-platcau region), ^ (the ^-plateau region), and 1 (the 
saturated magnetization region); the Bjj and Bld regions are, 
respectively, the 'Haldane-type- ^-plateau' and 'large- D-typc-^- 
plateau' regions. In other regions, m increases continuously with 
the increase of H/Ha. 

ues of H giving the ^-plateau, respectively, we extrapo- 
late the corresponding results for A'' = 36, 64, and 84 to 
N^oo, assuming again the A/'-dependences of the finite- 
size values are quadratic functions of A''"^. The magne- 
tization curve thus obtained for the cases of J =—1.0 
and D = 3.0 and of J = —10.0 and D = 3.0, which are, 
respectively, in the 'Haldane-type-^-plateau' and 'large- 
-D-type-^-plateau' regions (see Fig. 2), are depicted in 
Fig. 3. 

The magnetization phase diagram on the H/Hs ver- 
sus J plane for a given value of D is obtained by plotting 
iJo, Hi, H2, and Hg as a function of J. The results for 
the cases of D = 0.0 and _D = 3.0 are shown in Fig. 4. Fig- 
ure 4(a) demonstrates that in the case of D = 0.0, the 
^-plateau, which is of the 'Haldane-type', appears when 
-1.69 ±0.01 <J< 0.0, and Fig. 4(b) shows that in the 
case of £' = 3.0, the i-platcau disappears at the criti- 
cal point Jc^''"^-' = —3.21 ± 0.01 for the 'Haldane-type- 
l-plateau'-'large-f-type-^-plateau' phase transition (see 
Fig. 2). 

4. Concluding Remarks 

We have investigated the finite-field ground state of the 
S = 1 antiferromagnetic-ferromagnetic bond-alternating 



chain described by the Hamiltonian H [see eqs. (1-3)] 
with J<0 and — oo<D<cxd, focusing our attention 
mainly upon the ^-plateau appearing in the magneti- 
zation curve. We have found that two kinds of magne- 
tization plateaux, i.e., the 'Haldane-type-5-plateau' and 
the 'large- D-type-^-plateau' [Fig. 1] appear. We have de- 
termined the 5-plateau phase diagram on the D versus 
J plane [Fig. 2], applying the TBCLS methods.^^'^ We 
have also calculated, by means of the DMRG method,^ 
the ground-state magnetization curves for a variety val- 
ues of J with D fixed at D = 0.0 and 3.0, two examples 
of which being presented [Fig. 3]. Using the results for 
the magnetization curves, we finally obtained the mag- 
netization phase diagrams on the H/Hs versus J plane 
for D = 0.0 and 3.0 [Fig. 4]. 

For the purpose of exploring the effect of frustration, 
we axe now studying the case where the antiferromag- 
netic next-nearest-neighbor interaction term J2 J2e ■ 
(J2>0) is added to the Hamiltonian of eqs. (1-3), 
assuming either J > or J < 0. It is noted that when 
J = 0, the system is reduced to the S —1 antiferromag- 
netic ladder, where the ratio of the leg interaction con- 
stant to the rung one is J2. The results will be published 
in the near future. 
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